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Abstract 


A method is presented for calculating the aerody- 
namic heating and shear stresses at the wall for tangent . 
ogive noses that are slender enough to maintain an at- 
tached nose shock through that portion of flight during 
which heat transfer from the boundary layer to the wall 
is significant. The lower entropy of the attached nose 
shock combined with the inclusion of the streamwise ... 
[pressure gradient yields a reasonable estimate of the 
actual flow conditions . Both laminar and turbulent 
boundary layers are examined and an approximation of 
jthe effects of (up to) moderate angles-of-attack is in- 
cluded in the analysis. The analytical method has been 
programmed in Fortran IV for an IBM 360/91 computer. 

Symbols < 

C f = friction coefficient local (-) 

C = specific heat of air (Btu/lbm °K) 


(dv/dx) 


( 0 ) 


= Newtonian velocity gradient at the stagna- 
nation point on a circular nose (sec -1 ) 
(see Eq. 19) 


f" 

_ .( w ) 


= velocity gradient parameter from Refer- 
ence (6) (see Eqs. 27 and 28) 

g; — acceleration . of gravity (32.174. £t/seo2) 
h = enthalpy (Btu/lbm) / 

, h* = reference enthalpy (Btu/lbm) \ • -p.f r . 
’H* 0) = defined by Equation (18) 

H* = defined by Equation (17) 

C x ) . 

k ■■= coefficient of thermal conductivity • :/ 
(Btu/ft sec °K) . ., v 

N = a constant; set = 0 for one-dimensional 
flow; set = 1 for axisymmetric flow (= 1 
throughout this program) 

Nu = Nusselt Number (-) 

p (i) = pressure at stations 1 through 15.(atmo- . 
spheres, except where otherwise noted) 

\ 

P r = Prandtl Number (-) 

PXPO(l) = the local to stagnation point pressure 

ratio at each of the 15 station locations (-) 

A 

q = heat transfer rate (Btu/ft 2 sec) 

QRATL = ratio of heat transfer with cross flow to that 
without cross flow for a laminar boundary 
^ X layer 7 (see Eq. 29) (-) 


R e ° r R e 


QRATT = ratio of heat transfer, with cross flow to that 
without cross flow for a turbulent boundary 
layer (see Eq. 30) (-) 

R = tangent ogive radius of curvature (see Fig- 
ures 1 and 2) (ft) 

= local Reynolds number (Eq. 25) (-) " , 

= local momentum thickness Reynolds number] 
(Eq. 26) (-) ' 


(*> 


R„ 


-( 0 ) 


Rj^ = spherical nose radius (for calculations of 
1 blunt body 1 stagnation point heat transfer 
rate) (ft) 

r(i)or r (X) = the flow deflection distance defined by Eq. 

(11) and shown in Figures 1 and 2 (ft) 

V = velocity (ft/sec) 

X = tangent ogive longitudinal dimension (shown 
in Figures 1 and 2) (ft) 

X (I) = surface coordinate distance along stream- 
line from nose tip to each station (ft) 

j 

Y = tangent ogive base radius, shown in Figures' 

1 and 2 (ft) j 

I 

’ a = ALPHA = the vehicle angle- of- attack (deg. 
or rad.) j ■ 

S = the flow deflP.C.t.inn pngl© (oonioal flow) Q-t 

which the nose shock becomes detached for j 
a given free stream Mach number (deg. or 
• rad.) . . '■ . • • ; 

*/ 9 = the local surface deflection angle; also, the 
central angle turned by the tangent ogive 
radius, R, to define the complete ogive (see 
Figures 1 and 2) (deg. or rad.) : ’ 

6 = cone half-angle (deg. or rad.) 

= viscosity coefficient (lbf sec/ft 2 ) 


P 

P 




<£(!) 


= density of air (Ibm sec 2 /ft 4 = slugs/ft 3 ) 


= boundary layer shear stress at the wall / 
.• (lbf /ft 2 ) 7 ; ..j. 

= angle defined by Eq. (10) and shown in 
Figure 2 (deg. or rad.) 

Subscripts 


r-\ 


e = local, extemal-to- the- boundary layer value 
Lam = considers a laminar boundary layer 


rec = evaluated at recovery conditions 





-/ 


^ r e f = evaluated at reference conditions (see . ] 

superscript) ■'*' • 1 

Turb = considers a turbulent boundary layer 

x = at a position X feet from nose tip along a 
surface streamline (same as e) 

w = evaluated at local pressure and wall 
temperature 

o = at stagnation point for a spherical nose of 
Radius, R M 

1 N 

co = free stream!(ahead of nose shock) value 
Superscript 

* = property evaluated at local pressure and 
reference enthalpy 



Figure 1. Geometry of Tangent Ogive Nose 


Introduction 


The generally favorable aerodynamic characteristics 
of the tangent ogive in supersonic and hypersonic flow re- 
sult in the common use of this configuration for sounding 
rocket noses. Accordingly, an analytical method for cal- 
culating the aerodynamic heating on such configurations 
has been devised, combining basic analytical methods 
which are well known with some which are less common 
and with certain basic assumptions. These methods, 
while approximate in nature, yield results which have 
proved to be adequate for the design of both the structure 
of the rocket nose and the protection of payload items 
within. The entire analysis described here has been 
programmed in Fortran IV for an IBM 360/91 system 
(Reference 1). 

The slendetness of the ogives of interest results in 
an attached nose shock wave through periods of super- 
sonic and hypersonic flight during which significant aero- 
dynamic heating is experienced. The low entropy-jump 
across the oblique shock wave as opposed to the entropy-, 
jump across the normal shock wave associated with "blunt 
bodies" results in an increase of heat transfer to the , 
ogive for a constant flight condition. This is similar to 
the case of the cone heating as compared to that on a 
blunt, axisymmetric body. However, unlike the cone case, 
the ogive body has a definite (first order) pressure 
gradient along the surface streamlines. A blunt body 
analysis is treated in Reference (2) and a conical body 
analysis in Reference (3). The present analysis con- ; j-~ ; 
siders the in-between (tangent ogive) case in which the 
nose shock is oblique but there is a body pressure grad- 
ient. The effects of moderate angle- of- attack (local body 
angle plus angle-of- attack of 30 to 35 degrees) are 
approximated. 

Theory 

\ 

The theory is derived from a combination of the 
. analytical methods of References (2) and (3) with several 
new approximations and assumptions. The pertinent 
geometry along with the most important items of nomen- 
clature are shown in Figures 1 and 2. The tangent ogive 
and flow geometry are completely defined by the param- 
eters x , y and a (identified in Figure 1) in con- 

ju^ction with the g eometric equations (7 through l2)^The 


effects of (up to) moderate angles-of-attack are accounted 
for by assuming the local flow to be similar to that on a 
cone of half-angle (equal to the ogive local surface angle) 
at angle- of- attack. The applicable free stream conditions 
are derived either from the vehicle altitude and an ap- 
propriate ARDC atmosphere or by electing to define the 
jpr e-nose shock air properties by specifying two thermo- 
dynamic variables - the pressure and the temperature — 
jthen obtaining all the other properties from the real gas 
[(equilibrium) Mollier approximations of Hansen (Reference 
|4). !■ 


For simplicity, a fixed number of body locations are j * 
Specified for each problem. Either of two procedures can 
jbe adopted for defining the local flow conditions at these 
Specified 15 body points. In the first case, the local presj j 
|sures at the calculation points are defined by Newtonian ! 
'approximations or from experiment (if available). The \ 
'entropy behind the attached nose shock (assumed to be ! 
jconical) is then calculated and the local, external-to- ; 
boundary layer properties are defined by isentropically 
expanding to the given local pressures. This assumes that 





ithe er.tropy is constant at the post-nose shock value over 
the entire body and results in an Over-prediction of the 
available energy, hence, heat transfer rate to the wall. 

The second method obtains the post-nose shock entropy 
and pressure (as in the first method) and taking these | 
data as initial values, uses a Prandtl- Meyer expansion | 
through A 6 degrees to get the pressure at point 2 (Figure 
1). This pressure and the point 1 entropy are then used 
with the air properties of Reference (4) to define all ex- 
ternal- to-bounday layer properties at point 2. The local 
surface angle at point 2 is then considered to be the cone 
half-angle and a new cone f external- to-boundary layer) 
entropy at point 2 is calculated and used in the same 
manner as just described to expand by an appropriate 
Prandtl- Meyer routine through A 6 degrees to obtain the 
point 3 pressure at the point 2 entropy. In this way, the 
properties at each point on the ogive are approximated by 
using the entropy of the immediate upstream point. 

'• This method represents an attempt to approximate 
the entropy gradient across the shock layer. Note that 
the entropies so derived are conservatively lower the 
farther aft one goes on the ogive. Moreover, the greatest 
conservatism in the predictions of the local heat transfer 
rate occurs at point 11. For all points downstream of 11, 
the entropy is left at the point 11 calculated value and the 
pressures are assumed to be the arithmetic average of the 
free stream (ahead of the nose shock) and the point just 
upstream of the point being studied. Mathematically: 


SR (11) = SR (12) = SR (13) = SR (14) = SR (15) (1) 


| Geometric Calculations ; 

i 

! The programming of the analysis is greatly simplified 
by the fact that a fixed number of stations on a fixed geo- • 
metry (the tangent ogive followed by a cylinder) are j 
treated. In specifying the magnitude of x max and y max and 
jthe angle-of-attack, a , all necessary geometric input has 
been supplied and one is able to generate the geometric j 
details required by means of the following equations (the j 
nomenclature of which is defined in Figure 2): Given j 
X max , y max and a , j 


x 2 .+ y 2 


6 - Tan 1 


-^4 = Tan'* fey-, : 

L R -yJ v {_x 2 — y 2 . 


A0is defined as 


a a = o.i ( 0 ).. ~ ■ J .■ • .. . (® 

In order to calculate the "flow deflection" distance, r(i), 
define 

: ' , Z> 

4> (i) = PHI (i) = 6 + a (K 


P (11) +P (10))/2- 
P (12) =(P ro +P (ll))/2. 


and thence (from the right triangle of Figure 2) 


r (i) = 2 R sin sin 0 (i) (H)j 


P (13) =f> 0D +P (12))/2. 
P (14) =(p a + P(13f2. 


P (15) =P a 


The surface coordinate distance (boundary layer build-up j 
distance) is 


.X(i) =jR 




Of course, the Prandtl- Meyer expansion is used onlyj 
up to and including point 11. From point 11 to points 12, j 
13, 14, and 15, the entropy is constant and isentropic ex- j 
ipansion to the local pressures indicated in Equations (2) 
■through (6) defines the local flow properties. The calcu- 
lated pressure is not allowed to go below 80 percent of 
ithe free stream pressure - as would sometimes result 
from the Prandtl- Meyer expansion technique described. 

It should be emphasized that while the entropy values 
for the second method are low (yielding conservatively 
high heat transfer rates), the pressures that result from 
the Prandtl- Meyer expansion technique are also normally! 
on the low side. This tends to decrease the predicted 
heat transfer rates. The two effects tend to offset each j 
'other. : 


Assumptions j 

The following basic assumptions are applicable to the 
analytical methods: j 

! ■ 

1 . The shock layer entropy gradient effect on the 
exte mal- to- b oundary layer flow properties around 
the ogive is approximated by assuming the local 
external flow to have originated just downstream 
of a shock wave generated by a cone of the same 
half- angle as the local surface angle (9 + a). This 
tends to predict increasingly lower (than the true , 
local value) entropies as one considers points I 

farther downstream on the ogive. The predicted j 

heat transfer, therefore, is expected to become 
i increasingly conservative as the farther down- 

stream stations are treated. 

: 


2. -In the alternate method in which, pressures are in- 

put and the entropy, is cohstant at the post-nose 
shock (point 1) value, the calculated entropy is ex- 
pected to be higher than the actual values for 
downstream points, again the error growing with 
distance downstream. Accordingly, lower heat 
transfer rates downstream should result from the 
entropy effect of the pressure-input option. Note 
that one cannot conclude from this that the heat 
rate distributions from this method are actually 
conservative because their magnitudes in this 
option are also highly dependent upon the magni- 
tudes of the pressures that are input. The above 
remarks refer to the entropy effect alone (as 
though the pressures by either method were equal). 

3. It is assumed that the increase in heat transfer 
rate caused by the thinning of the boundary layer 
(resulting from the cross flow) at any station on 
the ogive nose can be approximated by the use of 
Equations (29) and (30), substituting the local sur- 
face angle (0) for the cone half- angle, 6 c . The 
validity of this local similarity assumption has not 
been established independently at this writing. 


N = 0 for two-dimensional and N = 1 for a.xisymmetric 
flow (hence N = 1 for the case considered) and 


!H : 


( 0 ) 


2 P 


P 


(x) ( d V 
d x 


( 0 ) 


L V » ^ A (x) // ^’(0)lJ 


1/2 


(1 + N) 1/2 


(18) 


is the stagnation point value. 


Note that (dv/dx) (Q) is given the Newtonian (circular 
nose) value 


-- 1 


'dV\ _ /2 
dx/ 0 \ 


P (0) P (”> 


P, 


(0) J 


1/2 


(19) 


■The turbulent heating equation is (Reference 6): 


Surb^ 


(a; x) ) 7/1s c c ; (0) ) 2/3 x - 2 


Heat Transfer Rate and Shear Stress on Body 


i i 


Using Reynolds analogy, the laminar and turbulent shear 
stress and friction coefficients are calculated from: 


The boundary layer heat transfer rate and shear 
stresses at the wall are calculated by means of (a) the 
Eckert and Tewfik adaptation of Lee's momentum integral 
jequation (Reference 5) and the use of Reynolds analogy for ~ 
the laminar case, and (b) the Flat Plate Reference Enthalpy 
IMethod described in Reference (6) (also applying Reynolds 
analogy) for the turbulent boundary layer case. The pro- _ 
grammed equations are: ! 

I ... | 

I 1 . Laminar heat rate (ratioed to the spherical nose | 
stagnation point heat rate) (Reference 5) j 


^ ^Lam (P r( x) ) 


2/3 


q ( x ) _ 

q (0) 


k* ^ 

(X) 


h?.; 


h 

r ec 

- h <w> 


CP w 

(0) 

k* 

(0) 


_ H (0)_ 


h ( o') 

- h o’) 


Cp w 

(x) 


(13) 


[where 'the starred quantities are evaluated at local pres- 
jsure and reference enthalpy, ,h*. 


- h* = h ref = ^- +h(W) + P.22 ^ (h (0) - h (w) ) (14 


Vec = h (e) [1 - ( p ;> 1/2 + h (Q) C0 1/2 ] (L a minarB . L .) < 15 j 


LBm P(x) V (x) ( h rec ~ h (w>) 6 

! c _ ^ Oturb (P r( x) ) 2/3 

fTurb ^(x) V (x) ( h rec " h (w)) B 

! 

i 

r T La m = °- 5 C f Lam ^(x) V (,) 

r- , =0.5C. p, v 2 

Turb f Tu rb (*) ( x ) 

The local Reynolds number is calculated from 


( 21 ) 

i 

( 22 ) 

(23) 

(24) 


:or 


Vec h (e) ^ ( p r ) 1/3 + h (0 ) (P r ) 1/3 ^ (TurbulentB.L./ 16 j 


R 


-'(*) V (x) X 


'(*) 


(25 


P, 


(x) 


'and the momentum thickness Reynolds number is (see 
Reference 2) j 


jin Equation (13), H*(x) is given by: 


i H oo 


(X) 


'<*> 


'( 0 ) 


(r(x)) N 


— 


~ p l *) 




f( 0) 

M (0). 


(^) C D LL • 

U f(w) ^(X) ^(x) 


(26) 


*(x) 


(c°) 


(r(x)) 2N d x 


(17) 


in which the Cohen and Reshotko's (Reference 7) velocity ! 
gradient parameter, f” , is curve fitted by the follow- j 
ing two equations (valid for favorable pressure gradients): 

‘ f V 

( a ) C f R e /N u > 2 j 

(27) 

f" w) = 0.0508 [C f R e /Nj 2 + 0.1332 [C f R e /N u ] ; 


1/2 
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TEMPERATURE 
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LOCATION 1 


OGI'Te-CYLINDER 1 
SHOULDER ' 1 
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, 4 ; Reference 2) is valid only "down to Iqcad Mach numbers ; 
somewhere between 2 and 3 . As with the cone shock cal- 
culations, the real gas Prandtl- Meyer analysis is backed 
, up, by a perfect gas routine (Appendix D, Reference 3) which 
" is used when failure of the real gas method occurs. Note J 
that in either the shock or expansion analysis for low 
^supersonic flow, the perfect gas assumption is quite valid. 


{Figure 3A. Sketch of the vehicle Nose for NASA TND889 | 


OGIVE- CYLINDER 
SHOULDER 


Cross Flow Corrections 


0.062 IN 

ST. STEEL WALL 


TC#1 
(inner wall face) 


TC#9 

(inner wall face) 


The derivation of the effects of cross flow on a cone j 
! at moderate angle- of- attack is taken from the work of j 
jj. Sternberg and is given in some detail in Reference (3). j 
The resulting equations are: 



QRATL = 


Figure 3B. Sketch of Black Brant VC Nose Showing 
Locations of Thermocouples 


*with crossflow 


^without crossflow I 

— 1 Laminar 


and 


r 1 

2 Tan (6 + a) 


2 -. 

Tan 6 

C 

\- 

L. - 

3 


1/2 


( 29 ), 


!(b) 0 < C f R/N < 2 


(28) 


QRATT = 




with cros s f low 


‘without crossflow 


- 1 Turbulent 


(30) 


.1 . .1 

— = 0.011627 [C f R e /Nj 2 1 + 0.25644 [C { R e /N u l - 0 . 08978 r 

i j 

1 = 0.85 (1 + k 2 )- 2 

i 1 ! 

i ' 

{where 

• 

Shock and Expansion Routines , 

1 

k 2 = 1.25 

Tan (& c + a) 


Tan 


l me real gas — - —j £ 

flow and the method of accounting for the pressure and 
'entropy gradients across the shock layer for conical flow 
are given in Appendix B of Reference (3). This hyper- 
sonic analysis (real gas in equilibrium) becomes invalid 
in the medium to low supersonic range so an alternate 
'supersonic conical shock calculation method (perfect gas 
l see Appendix C of Reference 3) is available upon failure 
of the hypersonic analysis. 


Similarly, a real gas (applicable for hypersonic Macti 
number ranges) Prandtl- Meyer expansion routine (Figure 



These equations are used to estimate the effects of cross ■ 
flow on the ogive nose. Both equations are solved at each 
station on the ogive assuming the local surface angle, 0 , ; 
to be the "cone half-angle," & c . Note that when the angle- 
of- attack, a, goes to zero, both QRATL and QRATT go to 
lone . F inally , it must be remembered that ev en when prop- 
er ly used (for cone flow) Equations (29) and (30) are validj 
only for (6 C + a) values up to approximately 35 degrees. J 
Sas the deflection angle plus angle-of-attack exceeds this | 
range, the cross flow begins to dominate the aeroheating j 
phenomenon and a better analytical prediction derives from 
a two-dimensional, blunt body method like that of Refer- j 
fence (2). Obviously, the crossover point for the applica- j 
bility of either theory is not clearly defined. 


Hemi-Spherical Nose Stagnation Heating j 

There are two reasons for including the calculation of 
the he mi- spherical nose calculation in this analysis. First, 

1 < • 4-"U m 


21 23 25 27 29 

TIME FROM LAUNCH -SEC 


Figure 4. Velocity and Altitude Histories for the Flight 
of the Vehicle of Figure 3A through the Significant 

Aeroheating Portion of Flight _ _ 


LUC UC11U ~ — 

■the post-normal shock flow properties axe required for the 
l aminar boundary layer heat rate calculation of Equation 
|( 13 ) and, in the same equation, the stagnation heat rate is 
'specifically required to redimens ionalize the heat ratio to 
get the absolute value of the local laminar heat rate. The, 
second reason lies in the fact that the stagnation point j 
heat transfer is often desired as a general heat transfer , 
parameter for evaluating the effects of trajectory param- 
eters, vehicle weight, etc., upon the vehicle thermal 
environment. / 

>5 


| I'Jote that the nose radius selected for the calculation 
Renters into the solution of Equation (13) only via the H* 0) i 
term (Equation 18) in the form of the stagnation point j 
[velocity gradient (Equation 19). The same nose radius is! 
also used in the stagnation point heating rate (Equation 31) 
so it is clear that all data relative to body points 1 through 
15 are totally independent of the value of RN. For this 
reason it is suggested that a value of RN = 1 foot be as- 
sumed because this is commonly used as a heat indicator. 

j 

The stagnation point heat rate can be calculated by the 
equation of Fay and Riddell (Reference 8) with the Lewis j 
number assumed equal to unity: 

i 
; 
l 

I 

(31) 


^stag. point = 0-76 g (P r ) ' 6 (P ( .) W 1 (V)'W 4 < h C 0 ) 


. . /dvV 7 

( w > \dxj 0 


Detachment of Nose Shock Wave 

The method of estimating the point 1 conditions by 
assuming the conical nose shock wave that would result 
from a cone of half-angle equal to the angle shown as 
(6+ a) in Figure 2 does introduce a possible problem. 

If (6 + a) at any given free stream Mach number is 
greater than some critical flow deflection angle, S cr , the} 
nose shock will detach and the analytical methods will 
become invalid. In order to protect against this possibility 
(which can arise from either too blunt an ogive, too large j 
an angle-of-attack, or a combination of the two), the 
critical conical flow deflection angle as a function of free! 
stream Mach number is taken from chart 5 of Reference ! 

I ! 

1(9). Thus, for each problem, the nose apex flow deflec- j_._ 
!tion angle (including angle-of-attack) and the free stream 
Mach number are known. The analyst uses the free 
stream Mach number in the appropriate region of Figure ! 

5 (Reference 9) to obtain the maximum allowable flow 
deflection angle, § cr . This angle is then compared with 
•the actual deflection angle ( 6 + a .) at t he nose and if. § cr j... 
-< (& + a), it is clear that the shock is detached and the 
present analysis is inapplicable. 


Comparison of Theory With Flight Data 


- 


In order to test the theory presented here with flight 
data, two cases are considered: a four stage vehicle re- 
ported in Reference (11) and Black Brant VC, Flight 
21.006 UG (Reference 10). 



i Figure 5. Laminar and Turbulent Heat Transfer Rates, 

; Recovery Enthalpy, Local Reynolds Number and Mo- 

j. mentum Thickness Numbers at Point ’JA n 


I. Comparison With Flight Data From NASA TND 889 
! (Reference 11) 

| The vehicle nose of Figure 3A was flown on a four 
[stage vehicle with a temperature history measured on the 
inside of the 0.032 inch inconel wall at the position marked 
["point A" on the sketch. The velocity and altitude histories 
[of the test vehicle through the significant heating portion ! 
[of flight (during which valid temperature data were re- I 
[corded) are shown in Figure 4. Using these trajectory j 
[data and the nose configuration of Figure 3A, the la m inar j 
[and turbulent heat rate, recovery enthalpy, and the local j 
[Reynolds number data of Figure 5 were calculated by the j 
[methods of this paper (using the computer program NQIH 
DW019 from Reference 1). The Prandtl- Meyer expansion; 
method of obtaining local pressures was used. 

These data were then input to a 10-element, one- 
dimensional structural heating program (NQLDW112, 
Reference 12) to obtain the temperature histories shown j 
as solid lines inFigure 6. Note that the cases of a fully j 
turbulent boundary layer and of transition at local Reynolds 
numbers of 2.8, 5.0 and 10.0 million are presented. The [ 
[flight recorded temperature data are shown as circled j 
points. 

i i 

I '■ ] 

In general, the agreement between theory and flight [ 

data is quite good. The data indicate that transition from ! 
[turbulent to laminar flow probably occurred at a calculated 
local Reynolds number of seven to eight million. It is 
[emphasized that this local Reynolds number is defined as: 



Re 


(TO 




( x ) 


(32) 


where the sub x values are taken at the outer edge of the 
boundary layer at point A. Inasmuch as the manner of ap- 
proximating the effects of the entropy gradient through the 
shock layer results in more or less fictitious values of - . 
the local entropy at any given point, caution must be used 
in comparing thd transition Reynolds number as defined 
by Equation (32) with transition Reynolds numbers from ,, 
other-sour ces-which. do not- make the same local- entropy" \ 

, ' i 

value assumptions. 

III. Comparison With Flight Data F rom Black Brant V C j 
| Flight 21.006 GT (Reference 10) j 

i | 

The vehicle nose of Figure 3B was flown on Flight j 

21.006 GT, a Black Brant VC, and two thermocouples were 
located on the inner surface of the 0.062 inch stainless ‘ 
steel wall at the positions indicated as TC #1 and TC #9 ; 
in the sketch. The trajectory used in the calculations is . 
taken from radar data and the velocity and altitude his- [ 
[tories are shown in Figure 7. The digital program . 
[(NQLDW019) was used to derive the heat transfer rate ; 
data of Figures 8 and 9 for TC #1 and TC #9, respectively. 
[These data, again using the structural heating analysis j 
[of Reference (12), result in the temperature predictions I 
((solid lines) of Figures 10 and 11. 

i . 1 * * 

i 

j 

If transition from turbulent to laminar flow is as- 
sumed to occur at a local Reynolds of ten million at TC 
#1, the agreement between theory and flight is seen / 
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Figure 6. Comparison of Flight with Theoretical Data j 
Using the Theory of NQLDW019 to Predict the Heat j 
Rate History 
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(Figure 10) to be extremely good. Similarly* the assump- 
tion of a transition Reynolds number of 15 million at TC 
#9 results in good agreement between theory and 
measurement. 

Transition Reynolds Number 

It cer tainly comes as no surprise that the ability to 
predict the heat transfer rate is primarily dependent upon 
the ability to predict the nature of the local boundary 
layer. This holds true whether the analytic approach be 
approximate (as in the present case) or involves a com- 
plete numerical solution of the boundary layer equations. 
Figures 6, 10 and 11 indicate that the present analysis is 
adequate if the transition can be predicted. Accordingly, 
'an empirical method for predicting the transition Rey- j 
jnolds number is sought by the simple expedient of record- 
ling experi mental valu es as they are^ inferred by j 
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Figure 8. Black Brant VC Flight 21.006GT Heat 
Transfer Rate at TCI 
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Figure 9. Black Brant VC Flight 21.006GT Heat 
Transfer Rate at TC9 

comparison of the heating analysis predictions with the 
Qight temperature data as has just been demonstrated. 


Such a comparison is shown in Figure 12. Of course, 
with the scant amount of data available the ability to 
predict transition Reynolds numbers for aU tangent ogives 
at all flight conditions is almost totally absent. On the | 
pther hand, it should be borne in mind that the flight re- ; 
gimes, degrees of bluntness of the ogives and general j 
body surface conditions for large families of sounding ! 
rockets are sufficiently restricted as to make such an 
'empirical approach quite practicable. In the present case 
[(Figure 12), only thin- wall, metallic pgives axe represented. 
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Figure 10. Black Brant VC Flight 21.006GT Tem- 
perature of 0.062" Stainless Steel Wall at TCI 



It is quite likely that low conductivity noses (i.e., fiber- 
glass phenolic) would have appreciably different ratios 
of external surface to local recovery temperatures and 
thence would be expected to require a separate curve 
such as that of Figure 12. Over a period of time, then, 
the use of the present theory should reasonably quickly 
allow an investigator to attain an adequate ability to pre- 
diet the effects of aerodynamic heating on any tangent 


ogive nose. 


Conclusions 

The following conclusions are . drawn from the work 
reported: 
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Figure 12. Summary of Transition Local Reynolds Num-| 
ber Which Best Matches Flight Temperatures with Pres- j 
ent Theoretical Predictions vs Ratio of Point Location on. 
| Centerline to Centerline of Ogive from Two Flights 


1 . The tangent ogive aeroheating analysis methods 
described are capable of providing adequate pre- 
dictions of the aerothermal effects on tangent 
ogive noses of typical sounding rockets. 


J 


2. The accuracy with which these predictions can be 1 
made is dependent primarily upon the accuracy j 
with which transition from turbulent to laminar i 
flow (in ascent) or laminar to turbulent flow (in 
re-entry) can be predicted. The extreme com- 
plexity of the boundary layer transition phenomenon 
suggests that the empirical approach is the most 
practical at this time, particularly when large ; 
numbers of similar vehicles are involved as is 
often the case with sounding rockets. 

i 

3. Only a relatively small amount of data for checking 

the validity of the conical cross flow analysis j 
(incorporated in the present method by the assump- 
tion of local similarity) has been obtained to date j 
so caution is advised in its use until further j 
corroboration is available. j 
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